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Abstract 

It is interesting to determine the fraction of the gluonium component in rj and r]' which has been 
under serious discussion for many years. Measurements on different decay and/or production modes 
were employed in literatures, thus larger uncertainties were unavoidable. In this paper we suggest 
to determine the mixing angles r] — rj' — G using the data of semileptonic decays of D and 
Dg. We extract the mixing angles (p' , (j)G and the model parameters simultaneously. Thanks to 
the new measurements carried out by CLEO Collaboration and there are sufficient decay modes 
to determine both the model parameters and mixing angles. The mixing angles from data are 
(f)' = (41.5 lb 2.0)° and sin c/ic = 0.00 it 0.36. Even though the central value of sin (j)G is still zero, an 
upper bound is set. Moreover, as suggested in literature, r7(1405) is a glueball candidate whereas in 
our picture, ?7(1405) may be identified as G with glueonium being its main content. Using all the 
model-parameters obtained above, we estimate the branching ratios of Df{D^) — )• Ge+i^e where G 
is identified as ?7(1405). 
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I. INTRODUCTION 



For a long time the mixing oi rj — rj' is of a great theoretical interest because it concerns 
many aspects of the underlying dynamic and the structure of pseudoscalar meson. One can 
investigate this mixing from two distinct schemes, namely they may be a mixture of the octet 
and singlet of the flavor SU(3), or ^[uu + dd\ and ss that are the mass eigenstates if we can 
assume an u-d degeneracy. The two schemes reflect different understandings of the essential 
physics. Moreover, the general principle of QCD implies that a gluonic degree of freedom 
of ^ may be involved in the physical states r] and ?7'[l-ll|. The QCD anomaly indeed 
may induce such a mixing. It would be crucially interesting to make a definite judgement 
if there is a sizable gluonic component in rj and 77' from either an underlying theory or 
phenomenology. Obviously even though the QCD anomaly is well formulated, it is hard 
to find its contribution to physical eigenstates because the matrix elements of the anomaly 
operator are dominated by the non-perturbative QCD and not reliably calculable so far. 
Therefore, one should turn to the phenomenological studies from which one may gain much 
information about the non-perturbative QCD, i.e. the mysterious aspect of the successful 
theory. In this work, we are going to determine the fraction of the glueball component in rj 
and 77' along the line. 



Some researches 12|-|l4j were done on this topic. Generally the authors obtained the 



mixing from different decay or production processes 15|. It can be conjectured that different 
modes are governed by different physics mechanisms besides the mixing which is the goal of 
our investigation, the uncertainties are not controllable, therefore the errors would be large. 
Instead, we suggest to use the data of semileptonic decays of D and Ds where the reaction 
mechanism should be the unique. Thus one can expect that in such a way, uncertainty 
could be greatly reduced. In the theoretical calculations, we employ the light-front- quark 
model (LFQM) which has been widely used for calculating the hadronic matrix elements. In 
the model, there are two independent model parameters (see the text for details) which are 
related to non-perturbative QCD and should be fixed by fitting data. 

In o ur p revious work jl6|, with the data of three independent measurements on BTZ{D^ — > 



17| . ElZ{Ds —)■ rje'^Ve) and BlZ{Ds 'ri'e'^Ve)uBi we determined the r} — rj' mixing 



(39.9±2.6(ea;p) ±2.3(t/ie))° and predicted the branching ratio Bn{D+ rj'e+i^e) = (2.12 ± 



0.23{exp) ± 0.20(t/ie)) x 10-^(l6| which is consistent with the current data Bn^^D+ 



r]'e+Ue) = (2.16 ± 0.53 ± 0.07) x 10"^ measured by CLEO collaborator recently(l9|. It is 
noted that there are three free parameters if the gluonic degree of freedom in rj and r]' is not 
accounted, namely two are the model parameters and another is the mixing angle. Instead, if 
one needs to consider the extra mixing between quark states with gluonic degree of freedom, 
there are four free parameters and at least four independent measurements are necessary. 

Even though in the framework the assumption of null component of gluonium state is 
adopted and the predicted branching ratio of B1Z{D^ — )• r/'e^Ue) is consistent with data, one 
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still cannot confirm that there indeed is no gluonium component in r] and 77' because the 
experimental errors are relatively large. Considering the present data, one can conclude that 
the central value of the fraction of gluonium in 77 and r]' is consistent with zero, but while 
taking into account the error tolerance, instead, we would only able to obtain an upper 
bound of its fraction. Fortunately, thanks to the new measurements carried out by the 



CLEO collaboration [19|, the data for all the four decay modes are available, which enable 
us to simultaneously fit the two model parameters and the two mixing angles. Thus we may 
determine the fraction of the gluonium component in r] and r]' ( concretely the upper bound 
of the gluonium component). 

After this introduction, we describe the working framework and provide all necessary 
formulas. In the following section, we present our numerical results and the last section is 
devoted to the discussions and conclusion. 



II. THE MODEL CALCULATION 



In the theoretical framework where gluonium component is not accounted, the mixing 



matrix of 77 — 77' is set as 
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where rjg = -^{uu + dd) and r]s = ss and it is noted that both of them are not SU(3) singlet. 
We referred this mixing as scenario-I in this work. 

Since the QCD anomaly causes the mixing between 77 and r]', there is no any rule to forbid 
a mixing between the quark states and a glueball state of the quantum number 0""*". 

As one extends the picture to involve a gluonium component, a new scenario which we 
refer as the scenario-II, was suggested in Refs. as 



(v] 




^ cos 0' 


— sin 0' 


\ 








sin 0' cos 0G 


cos 0' cos 0G sin q 




Vs 


\g) 




^ — sin ^' sin 0g 


— cos 0' sin 0G cos ^ 


) 


\9 J 



(2) 



where \q) = | gluonium) is a pure gluonium state and the physical state G was identified as 

As discussed in the introduction, once there are data on — > rje^Ue available in addition 
to the other three modes, we will be able to extract (pc directly. Moreover new BTZ{D^ — j- 
r^e+z/e) = (11.4 ±0.9 ±0.4) x 10 
ratios of BlZ{Da — )• rj{rj')e^h', 



^ deviates from the old datum and the combined branching 



23| demands a redetermination of the mixing angles in 



Our strategy is that the mixing angles 



be and the model parameters 



'?('?') 



m 



the wave functions of t] and i]' are simultaneously determined by fitting solely one type of 
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FIG. 1: Feynman diagram for meson transition amplitude 



data i.e. we theoretically calculate the branching ratios of rje^Vei v'^^'^e, 

Ds — )■ rje^Ve and Ds — ?■ rj'e^Ve and match them with data. Thus we "extract" the mixing 
angles directly from the semileptonic decays of D and Dg. 

We use the light front quark model (LFQM) to evaluate the hadronic transition matrix 
elements 24|-|34|. then obtain the decay widths of — )■ rj{rj')e^Vf. and Dg — ?■ ri{ri')e^Ve 
which are functions of 0', (pc /^^(,;') The transition diagram is shown in Fig. [H 

The hadronic matrix elements for P — )■ P transition can be parameterized as 



(3) 



where P represents the pseudoscalar meson, P'{P") is the momentum of initial (final) meson, 
q = P' - P" and V = P' - P". 

Functions f+{q^) can be calculated in the LFQM and their explicit expressions were 
presented as |25| . 
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where m[, m'l and m2 are the corresponding quark masses, M' and M" are the masses of 
the initial and final mesons respectively. All other notations can be found in the Appendix. 

Since the calculations are done in space-like region one need to analytically continue them 
to the time-like region. It is convenient to redefine the matrix elements as 

M'2 - M"2 



i,P{P")\A,\P{P')) = [V, -5 gJ Pi(?') 



q' 



-qMq'). 



(5) 
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The relations among the form factors are 



-.2 



Fi{q') = Mq'), Foiq') = Uiq') + (6) 
The form factors in the space-hke region are parameterized in a three-parameter form as 



l-a(^)+6 



Af2 



where F represents the form factor Fi and Fq. The parameters a, b and F{0) are fixed by 
performing a three-parameter fit to the form factors in the space-hke region. We then use 
these parameters to determine the physical form factors in the time-like region. 

The input quark masses are directly taken from Ref. j25| as = 0.26 GeV, = 0.37 
GeV, m, = lA GeV. 

We first need to determine the model parameters for D and Dg. They can be fully 
determined by fitting their decay constants which are related to their total widths. For a 
pseudoscalar meson the decay constant can be evaluated 



■fp = / dx2(fp'±—^^4:{m[x2 + m2Xi). 



The decay constants are experimentally measured as = 0.221 MeV, = 0.27 MeV. 
With them the model parameters for /J^^ and are fixed to be = 0.499 GeV, f3D, = 0.592 
GeV [sBl- Since the total widths are better measured, the model parameters of the heavy 
mesons D and Ds are determined with higher accuracy. 

Once the model parameters at the initial side are fixed, we would turn to concern the 
decay products. At first look, there are six free parameters /3^, (3^,, (3^, /3^,, (j)' and (pc to 
be fixed. It seems that there are not enough equations to determine all these parameters. 



However, as discussed in Ref. [16| two relations = /3'^, and (3^ = (3^, hold, thus the number 



of unknowns reduces to four. Matching the theoretically calculated branching ratios with the 
data for the four decay modes, we obtain all the values of (j)', (pc, (^'^{n')^ ^riiri')- '^^^ mixing 
angles are fitted to be 0' = (41.5 ± 2.0)° and sin^c = 0.00 ± O.36(sin2 0G = 0.00 ~ 0.13) 
where the errors are from the experimental side. The model parameters are = 0.35 

GeV and = 0.59 GeV. 

With these parameters, assuming ?7(1405) to be the physical state G whose main content 
is glueonium g, we are able to estimate the branching ratios of decaying into eta(1405) 
via its qq components and we obtain B7l{Df — Ge'^z/g) = ~ 8.6 x 10~^ and BTZ{D^ — )■ 
Ge~^Ue) = ~ 1.1 X 10~^. Accurate measurements on these semi-leptonic decay modes may 
tell us if the main content of ?7(1405) is glueonium (i.e. a glueball candidate) and further 
test the Scenario II for the mixing. 
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III. DISCUSSION 



In Ref.[9| by fitting the data of several radiative decay modes such as w — > 777, p — 777 
and u — )■ 7r°7 sin^0G = 0.115 ± 0.036^ and cp' = (40.4 ± 0.6)° was fixed if one sets sin^^c a-s 
a free parameter. 

The value of (p' obtained from different kinds of experiments are consistent with each 
other in a reasonable error tolerance but the central value of sin^^c diverges over a range. 
The central value of sin^^c obtained from the semileptonic decays of D and Dg inclines to 
that the fraction of the gluonium state in rj and t]' is consistent with zero, but since there 
exists a relatively large uncertainty in the data, it is hard to conclude that it does not exist, 
i.e. sin^(^G = 0. It is urgent to improve the precision of measurement on these channels 
especially on D{Ds) r( ev^. 

In this work, we employ the LFQM model to evaluate the hadronic transition matrix 
elements. We fix the model parameters by fitting data, so that some theoretical uncertainties 
are involved in those parameters. Definitely the results are still model dependent because 
we employ a concrete model: the LFQM. But since all inputs are taken from the same 
source (i.e. the data of the semi-leptonic decays of D and Ds), one can expect that relative 
errors would be partly compensated, so that the model-dependence of the results is somehow 
alleviated. 

Because of absence of the final state interactions, the semileptonic decays have obvious 
advantages for determining the properties of the produced light hadrons, such as the structure 
of /o(980), 77 — 77' mixing and even a mixing of pseudoscalar mesons with glueball. 

Moreover, as suggested, assuming that the physical state ?7(1405) is the mixture G whose 
main content is glueonium, we calculate the branching ratios of BTZ{Df — t- Ge'^Ue) and 
and BTZ{D^ — )■ Ge^Vg). The semileptonic decay modes will be measured by the future 
experiments. Then the proposed scenario II for the mixing of qq, ss and glueoinum g will 
be tested. 
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the error comes from the method. 
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Appendix A 



Here we list some variables appearing in the context. The incoming (outgoing) meson in 
Fig. [T] has the momentum P'^"^ = p'^^"^ + P2 where p[^"^ and p2 are the momenta of the 
off-shell quark and antiquark and 

p'+ = XiP'+, pt = X2P'+, 

Pl± = XlP± + P±^ P2± = X2P± - P'_L, (Al) 

with Xi and p'j_ are internal variables and xi + X2 = 1- 
The variables Mq, Mq, hp and A''( are defined as 

^,2 _ P'l + ^? ^ P'l + 

° Xi X2 ' 

= Jm(^ - (m'l - m2)2. (A2) 



.;.(M»_M?)^^^.'. (A3) 



where 



with v' - ^ - ""'^^'^ 

WKn p^ - 2 2X2M^ • 

N{=Xi{M''' - M^^). (A5) 
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